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ABSTRACT 

Numerical simulations of structure formation have suggested that there 
exists a good correlation between the halo concentration c (or the characteristic 
density 5c) and the virial mass Myj^ for virialized dark halo described by 
the Navarro, Prenk & White (1995) density profile. In this Letter, wc present 
an observational determination of the c-Myj^ (or Sc-M^ir) relation in the mass 
range of ~ IO^^Mq < M^r <~ IO^^Mq using a sample of 63 X-ray luminous 
clusters. The best-fit power law relation, which is roughly independent of the 
values of and Ha, is c oc M^^'^ or 6c oc M~^'^, indicating n ~ —0.7 for a 
scale-free power spectrum of the primordial density fluctuations. We discuss the 
possible reasons for the conflict with the predictions by typical CDM models 
such as SCDM, LCDM and OCDM. 

Subject headings: cosmology: theory — galaxies: clusters: general — X-rays: 
galaxies 
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1. Introduction 

High-resolution simulations of structure formation have suggested that the virialized 
dark matter halos with masses spanning several orders of magnitude should follow a 
universal density profile (Navarro, Frenk & White 1995; NFW) 

where ps and are the characteristic density and length, respectively. The former is further 
related to the characteristic density 5c and the critical density pc = (3i/g/87rG')Z(z) of the 
Universe at redshift z through pg — ScPc, in which Z{z) = (1 + 2;)^(0m/^m(-s)), and flM is 
the cosmic density parameter. Although the NFW profile was first obtained in the SCDM 
model, subsequent numerical studies have shown that this density profile is independent of 
mass, initial density fluctuation or cosmology (e.g. Cole & Lacey 1996; Navarro, Frenk & 
White 1997; Eke, Navarro & Frenk 1998; Jing 1999; etc.). The two free parameters in NFW 
profile can be determined from the halo concentration c = r^ir/rg and the virial mass M^j^ 
if the overdensity of the dark matter with respect to the average value is Ac, 

" 3 ln(l + c) -c/(l + c)' ^ ' 

It has been well established based on numerical simulations that there exists a good 
correlation between the halo concentration c (or the characteristic density 5c) and the 
virial mass M^j^ for any particular cosmology (e.g. Navarro, Frenk & White 1996, 1997; 
Salvador-Sole, Solanes & Manrique 1998; etc.). For example, given Ac = 200, in the mass 
range of 3 x lO^^M© < M^ir < 3 x lO^^M©, the c-Mmr relation can be well fitted by a 
power-law function: c = a{M^ir/MQ)\ with (a, 6) = (891,-0.14) and (186,-0.10) for 
SCDM model {Qm = 1-0, = Q, h ^ 0.5 and ag = 0.65) and ACDM model (Om = 0.25, 
Qa = 0.75, h = 0.75 and as = 1.3), respectively. While the NFW profile and the c-M^ir 
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correlation are only determined from simulations for halo masses as small as M^ij. — 10^^ Mq 
or as large as M^r = IO^^Mq, there is no reason to beheve that these results would not be 
valid for halos which are one order of magnitude lower or higher in mass (Burkert & Silk 



On the other hand, observationally we can only measure the distribution of baryonic 
matter rather than that of dark halo. It is crucial to link numerical predictions with 
astronomical observations. In a fully virialized system, the distribution and motion of 
baryonic matter can be determined by the underlying gravitational potential of the dark 
halo, if the velocity dispersion or temperature profile can be well measured. It is thus 
possible to work out the dynamical properties of the dark matter halo such as M^r, 
Tg, r^jir, c or Sc from the observed distribution of luminous matter in a gravitationally 
bound system. This will eventually allow us to examine observationally whether there 
is a correlation between c (or Sc) and M^r. Finally, a comparison of the observationally 
determined and numerically simulated c-Mmr (or 6c-Mmr) correlations may provide a useful 
way to distinguish different cosmological models. In this Letter, wc will make an attempt 
for the first time to derive the c-My^r and 5c~^vir correlations on cluster scales by fitting 
the observed X-ray surface brightness profiles of clusters to those predicted by the NFW 
profile as the cluster dark halos under the assumption of isothermality. Note that there 
is a striking similarity between the distribution of intracluster gas tracing the dark halo 
of the NFW potential and the conventional (3 model. Makino, Sasaki & Suto (1998) have 
explicitly shown that the NFW profile via isothermal hydrostatic equilibrium results in an 
analytic form of gas distribution 



in which a = ATrGfinipPsr'^/kT, and /x = 0.585 denotes the mean molecular weight. Briefly, 
our task is to determine the best-fit parameters of a and for an ensemble of X-ray 



1999). 
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clusters in terms of eq.(4), and then derive the relevant parameters Ps and 5c in conjunction 
with the X-ray temperature T. Finally, we examine the c-M^i^ and Sc-M^ir correlations by 
solving eqs.(2) and (3). 



2. Sample 

Essentially, we select our cluster sample from two ROSAT PSPC cluster catalogs: the 
36 high-luminosity clusters by Ettori & Fabian (1999; EF) and the 45 nearby clusters by 
Mohr, Mathiesen & Evrard (1999; MME). The EF clusters have high X-ray luminosity 
Lx > 10"^^ erg s~^, and a great fraction of them also have relatively high redshift z ~ 0.1-0.3. 
The best-fit values of a and for all the 36 EF clusters have been given by EF for a 
cosmological model of Hq — 50 km s"^ Mpc~^ and Qm — A conversion of into the 
values in different cosmological models should be properly made, if needed. The MME 
clusters are taken from an X-ray flux limited sample, and thereby located at relatively 
small redshift. We perform the fit of the observed X-ray surface brightness profiles 
of the MME clusters to the theoretical prediction Sx (x J n^gagM according to thermal 
bremsstrahlung, where the integral is made along the line of sight. A cross-identification 
reveals that 16 clusters are listed in both samples, for which our fitted values of a and are 
fairly consistent with those obtained by EF. We further require that the X-ray temperature 
should be obscrvationally available, and we take the temperature data from Wu, Xue & 
Fang (1999; and references therein). The finally merged EF and MME sample to be used in 
our following analysis contains 63 clusters. 



- 6 - 



3. The c-M^ir and Sc-M^ij. correlations 

We confine ourselves to a fiat cosmological model with VLm + = 1- In this case, 
the density contrast depends on the value of VLm and can be approximated by (Eke et 
al. 1998) Ac = 178^]M(-^)°•^^ and Z{z) = (1 + zf{l + zVLm + [(1 + z)-"^ - 1]VL^} and 
^m{,z) = Qm{'^ + z)^/Z{z). For a given cosmological model {Qm, ^a), we convert our best-fit 
values of a and into 6c and then obtain the concentration parameter c and the virial 
mass Myjr by numerically solving eqs.(2) and (3). We display in Fig.l an example of the 
resultant c and 6c versus M^i^ for Qm = 0.3 and = 0.7. It is apparent that there exists a 
strong correlation between c or 6c and M^r- We employ the Monte-Carlo simulations and 
the fitting method to obtain the best-fit c-Mmr and 6c-Mmr relations which are assumed 
to be a power-law. This enables us to include the measurement uncertainties in both axes. 
Note that for the EF clusters we have not accounted for the uncertainties arising from the 
fitted parameters of a and since EF provided no information about their error estimates. 
The results for a set of cosmological models are listed in Table 1. Additionally, we have 
tried the orthogonal distance regression technique ODRPACK (e.g. Feigelson & Babu 1992) 
and reached a steeper power index. For example, in the case of VLm = 0.3 and I^a = 0.7 we 
find c = lOi^-47±i-O3(M,,,/M0)-o-92±o.o7 ^ W^-^^^^-''\M,„/Mq)-^-'^''^'^-^'> (see Fig.l). 

Although ODRPACK can simultaneously account for the scatters around M^j^ and c (or 
6c)-, the goodness of the fit has not been improved for our problem in the sense that the 
reduced for both c-Mmr and 6c-M.mr relations are significantly larger than those obtained 
using the fitting method along with Monte-Carlo simulations. In Table 1 we have also 
given the power index for a scale-free power spectrum of initial density fiuctuations implied 
by our best-fit 6c-M^„ relation: 6c oc M'^^^^^"^ (NEW), which yields n ^ -0.7. 



EDITOR: PLACE FIGURE |T| HERE. 
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4. Discussion and conclusions 

In a virialized system, the distribution and motion of galaxies and intracluster gas (if 
their self-gravity is negligible) are determined by the underlying gravitational potential of 
the dark matter halo, provided that the velocity dispersion and temperature profiles are 
well measured. Therefore, one can probe the dynamical properties of the dark matter halo, 
though it is invisible, by using optical/X-ray observations coupled with the hydrostatic 
equilibrium hypothesis. In this Letter, we have made an attempt to derive the halo 
concentration c, the characteristic density Sc and the virial mass Myj^ for galaxy clusters 
characterized by the NFW profile from the observed distribution and temperature of X-ray 
emitting gas. Although the correlation between c (or Sc) and M^r has been well predicted 
from numerical simulations, this is the first time to determine the c-Myir and Sc-Myir 
correlations making use of the real data from astronomical observations. 

The correlation between c (or Sc) and Myj,. established in this Letter is applicable 
to massive halos in the mass range of ~ 10^'^ Mq < My^j. <~ IO^^Mq. However, we 
notice that the resultant slope —0.5) of the logc-logMyj^ relationship is significantly 
steeper while the spectrum (n ~ —0.7 ± 0.3) of the primordial density fluctuations is 
much flatter than the values predicted by typical CDM spectra for the mass halos with 
3 X lO^^M© < Myir < 3 X lO^^M© (Burkert & Silk 1999), especially on cluster scales where 
n Ri -2 (e.g. Henry & Arnaud 1991; Mathiesen & Evrard 1998; Donahue & Voit 1999; 
Mahdavi 1999; etc.). If our results are not a statistical fluke (Note the large dispersion 
of the X-ray data points in Fig.l), the above conflict may imply that we should abandon 
at least one of our working hypotheses: (l)Intracluster gas is in hydrostatic equilibrium; 
(2) Intracluster gas has isothermal temperature profile; (3)The self-gravity of intracluster 
gas is considerably small as compared with the contribution of the dark matter halo; 
(4)The NFW profile provides a precise description of the dark matter distribution. Yet, 



-8- 



further investigations should be made towards a robust constraint on the c-M^ij. and 
Sc-Mjjir relations before we come to a detailed study of the possible reasons for the reported 
discrepancy. 

We note from Table 1 that our best-fit c-M^ir and 6c-M,uir relations and the constraints 
on n are roughly independent of the cosmological models (f)^ and fl\). This property 
will be significant for distinguishing different cosmological models when combined with 
high-resolution numerical simulations. 

We gratefully acknowledge the valuable comments by an anonymous referee. This work 
was supported by the National Science Foundation of China, under Grant No. 19725311. 
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Table 1. The best-fit c-Mmr and 5c-Mmr relations 







c-M^ir correlation 






5c-Myir correlation 




n 


0.15 


0.85 


C=108-17±0.58(m„,,/Mq)- 


-0.50±0.04 


Sc 




-1.14±0.10 


-0.72 ±0.20 


0.30 


0.70 


c = lO**-23±O-«(M„i^/M0)- 


-0.51±0.04 


Sc 


= 1O2O-83±1-64(M„,,/M0)- 


-1.15±0.11 


-0.70 ±0.22 


0.50 


0.50 


C=lO8-29±O-69(M„i,/M0)- 


-0.51±0.06 


So 


= lO2O-89±l-76(M„i,/M0)- 


-1.16±0.12 


-0.68 ±0.24 


0.70 


0.30 


C=1O8-32±O-74(M,,,/M0)- 


-0.52±0.05 


5c 


= 1O2O-92±1-86(M„.^/M0)- 


■1.16±0.13 


-0.68 ±0.26 


1.00 


0.00 


C=lO8-35±O.81(M„i,/M0)- 


-0.53±0.06 


Sc 


= lO2O-93±1.99(M,i,/M0)- 


■1.17±0.14 


-0.66 ±0.28 
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Table 2. Cluster sample* 



cluster z T (keV) a (Mpc) 



A85 


0. 


,0559 


6 20+°-^° 


8 465+0-118 


n 360+0.021 

U.CSOU_Q 021 


A119 


0, 


,0443 


Vr I-Q+0.27 

0.0»_g 27 


12.30ll:r« 


2.585l°:355 


A262 


0, 


,0169 




6.71812:11^3 


094+°-°°* 


A401 


0. 


,0737 


8.00l°;fo 


a Sfi7+0.092 


n q3n+o-027 


A426 


0. 


,0179 


6 79+012 
"J- '»_0.12 


8.08410-108 


n i«i +0.012 

U-iOi_o.oi2 


A478 


0. 


,0881 


0.»U_g 35 


9 507+0-10'' 


419+0.018 


A496 


0. 


,0325 


4 13+008 


8 001+0-118 


195+0-014 


A520 


0. 


,2010 




11 35+0-0 


1 70+0-0 

1- "J-0.0 


A545 


0, 


,1530 


5.50l«-o 


12 52+0-0 


1 51+0-0 

i-oi-o.o 


A586 


0, 


1710 


a CI +1.15 


O 01 +0.0 

^-81-0.0 


26+°-° 


A644 


0, 


,0704 




o 7^1+0.158 
'1_0.158 


, -,„t.+0.051 
-^•-^ ' '^-0.051 


A665 


0. 


,1816 


8 26+0-90 


10.6910.0 


1 49+0.0 


A754 


0. 


,0535 


9 00+0-50 


1 n QQ+0.373 
lU.OO_g 373 


774+0.066 

' '*-0.066 


A780 


0. 


,0552 


o c:7+0.10 

■^■^'-o.io 


9 441+0.372 


n 9S5+0.046 


A1060 


0, 


,0126 


o 9^+0. 06 

■-■■^^-o.oe 


o c;.o+0.199 


314+0.018 
^••^1^-0.018 


A1068 


0, 


1390 


r 50+0.90 

'^■^'J-o.no 


9 72+0-0 


42+°-" 


A1367 


0, 


,0214 


3.5oi;^:« 


in as+O-696 

iU.«»_Q ggg 


1 Sfi-H+O.l™ 
1-«0.5_(,.179 


A1413 


0. 


,1427 




9.17+r 


0.57l°:° 


A1651 


0. 


,0825 


6 10+0-20 




0.5631°:°^^ 


A1656 


0. 


,0231 


o Qo+0.34 
»-3»_o.34 


22.3612-J67 


3 471+0.414 

•^■^'1-0.414 


A1689 


0, 


1810 


9.02l0:« 


10.93l«:l« 


715+°-°''4 

^-'l'^-0.074 


A1763 


0, 


,1870 


8-98li:84 


9.2010-0 


1 16+°-° 


A1795 


0. 


,0631 




10.0410-148 


0.462l°:°f5 


A1835 


0. 


,2523 


9 80+l-*0 


10.2210.0 


0.32l°:° 


A2029 


0. 


,0765 


8 47+0-41 


Q iqQ+0.133 
a.j.ao_o j^33 


n 3qn+o.o22 

U.J9U_g 022 


A2052 


0. 


,0348 




8.358l°;lS 


204+0 014 


A2063 


0. 


,0355 


3.68l°:ll 


8.19410:140 


3fi2+0-020 
U.ODZ_o 020 


A2142 


0, 


,0899 


Q 70+1-30 


q 44fi+0'357 
"-^^"-0.357 


636+0-091 


A2163 


0, 


,2030 


1 4 RQ+0-85 


9 16+0-0 


1 09+0-0 


A2199 


0. 


,0299 


4 10+008 


9 131+0 081 


oi,r+0.012 


A2204 


0. 


,1523 


9 20+1-50 


Q f;qQ+0.122 
O.OOO g -1^22 


0.236l°:°l« 


A2218 


0. 


,1710 


7 10+0-20 

'•^"-0.20 


10.321°:° 


0.9910:0 


A2219 


0. 


,2280 


12 40+0--''0 

1^-^^-0.50 


11.511°:° 


1 59+0-0 


A2244 


0, 


,0970 


8 47+0-43 

°-^'-0.42 


8.033l°:223 


35fi+0.037 


A2255 


0, 


,0808 


7 30+1-10 


24.801^:635 


5 991+1-982 
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Table 2 — Continued 



cluster z T (keV) a Va (Mpc) 





0. 


.0581 


7 CI +0.19 


1 Q OI +0.655 

10.^1_g ggg 


2, 


f,9f.+0.180 
^^"-0.180 


A2319 





.0559 


''■'■■^-0.15 


8 306+°-"^ 


0. 


on, +0.041 

"'-'-'--0.041 




0. 


2279 


11 10+1™' 


9 25+"-" 


0. 


64+"-" 


A2507 


0. 


.1960 




12 53+"-" 


2, 


63+"-" 


A2597 


0. 


.0852 




fi 9'iR+''-783 
D.ZOD_Q 


0. 


.(,.+0.129 
■^"^-0.129 


A2744 


0. 


.3080 


11 00+?„^° 


15 21+2-° 

J.U..^±_g g 


2. 


76+°-2 

■ '"-0.0 


A3112 





.0746 


4 24+°-^* 




0. 


90Q+0010 
^"»_0.010 


A3158 


0, 


.0575 




10 28+?-2S2 

J.U..^O_g 3g3 


1, 


093+9-081 


A3266 


0, 


,0594 


8 00+°-^° 


14 1 1+0.889 


2, 


.00+0.243 
■^■^"-0.243 


A3391 


0. 


.0553 


5 40+0.60 


y 9^^+0-297 
/.ZUO_g 


0. 


C-7Q+0.065 
■"-'"'-0.066 


A3^2fi 


0. 


.0114 


3 68+0-0'3 


6 985+''-°''4 

D.aOO_g g94 


0. 


|-,c7+0.004 
■^^'-0.004 




0. 


.0559 


^•*^-1.30 


9 146+0-380 
^•-'■^"-o.aso 


0. 


Q04+O.O79 
■^"^^-0.079 




0. 


.0475 


t- 19+0.20 
''•^^-0.20 


7 Q7c;+0.206 
'■o'0_0.205 


0. 


C7Q+0.040 
■C'9-0.040 


A3562 


0, 


.0478 


3.80+°-80 

— U.OU 


6.938+°°98 
— u.uyy 


0. 


070+0.022 
■ O 1 Q Q22 


A3571 


0. 


.0396 


6 7S+0-17 
D-'-^-O.lT 


Q 9C;n + 0.205 
205 


0, 


ggg + 0.040 

■"'-"'-0.040 


A3667 


0, 


,0566 


7 Q+0.6 
'■'-'-0.6 


o co-l +0.158 
«-°*^-0.158 


1, 


.yc+O.OSl 
■l'"-0.061 


A4038 


0, 


,0302 


o on + 1.60 

O.OU_g go 


c- 9+0.082 
O-84^_0.082 


0, 


047+0-004 

■"^' -0.004 


A4059 


0, 


,0478 


o Q7+0.12 

"^•^'-0.12 


8 645+0-l''8 

O.D40_g J9g 


0, 


009+0.023 
■ '^'^^-0.023 


AWM7 


0, 


.0176 


o 7c:+0.09 


8 756+0-19'* 


0, 


ot.7+0.024 
■•^^'-0.024 


Cygnus-A 


0, 


.0570 




7-i5iiro^; 


0, 


1f,1 +0.008 
■^"^-0.008 


IRAS 09104 


0, 


.4420 


8.501-4O 


10.09tH 


0, 


18+0° 
■i°-o.o 


MKW3s 


0, 


,0434 


3 00+"--'" 


8.422i;-- 


0, 


210-1-0.012 

■^-'■-'-0.012 


MS1358 


0, 


,3283 


7 50+4--'" 


14.29in 


1, 


■+°-o.o 


MS2137 


0, 


,3130 


/I O7-+0.38 
^■■^1-0.72 


11.48iro 


0, 


,0+0.0 

■-^^-0.0 


Ophia-A 


0, 


,0280 


9-ioiS;fo 


9.00410-JS8 


0, 


t.70+0.093 

■■-^'"-0.093 


PKS 0745-19 


0, 


,1028 


8.70ll:«° 


o qon+0 096 

0.9CSU_g g9g 


0, 


97(,+0.014 
■"""-0.014 


Tria-A 


0, 


,0510 


10 o^+o-s'' 


Q 977+0.088 
' -0.088 


1, 


042+0-023 
■"^^-0.023 


ZW3146 


0, 


,2906 


D.OO_g 34 


10.261°:° 


0, 


19+0-0 



*This table is only available in electronic form. 



-13- 



Fig. 1. — The c-M^^j. and Sc-M^ir relations for 63 clusters in the case of Qm — 0.3 and 
Qa = 0.7. The MME and EF clusters are represented by the filled and open circles, 
respectively. Note that the error bars for the EF clusters have only accounted for the 
temperature uncertainties. The solid lines are the best fitted power law relations to the 
data sets, while the dotted hues represent the results using the ODRPACK fitting method. 



